Coulomb interaction effects on nonlinear optical response 
in C6o 5 C70, and higher fullerenes 
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q ■ Nonlinear optical properties in the fullerene C 60 and the extracted higher fullerenes - C 70 , C7 6 , 
O ' 

C78, and Cs4 - are theoretically investigated by using the exciton formalism and the sum-over- 
states method. We find that off-resonant third order susceptibilities of higher fullerenes are a 
few times larger than those of Cqq. The magnitude of nonlinearity increases as the optical gap 
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qq ■ decreases in higher fullerenes. The nonlinearity is nearly proportional to the fourth power of 
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the carbon number when the onsite Coulomb repulsion is It or At, t being the nearest neighbor 
hopping integral. This result, indicating important roles of Coulomb interactions, agrees with 
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quantum chemical calculations of higher fullerenes. 
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It has been found that Ceo thin films show large optical nonlinearities [1-4] which are at- 
tractive from the viewpoint of scientific interest as well as technological applications. The 
magnitude of the third-order nonlinear susceptibility, XthgI^) = X^H^'i °°i 00 ■• w )i f° r third- 
harmonic generation (THG) is of the order of 10~ 12 esu to 1CT 11 esu. This large response is 
comparable to the responses measured for polydiacetylenes. The optical spectra of C70 [4] and 
higher fullerenes (C76, C 78 , C 8 4, etc.) [5,6] have been obtained. In order to explain the results of 
several interesting experiments, theoretical investigations [7-11] have been performed. We have 
applied a tight binding model [7] to C 60 , and have analyzed the nonlinear optical properties of 
C 60 . Coulomb interaction effects on the absorption spectra and the optical nonlinearity of C 60 
have also been studied [10]. We have found that the linear absorption spectra of C 60 and C70 are 
well explained by the Frenkel exciton picture [11] except for the charge transfer exciton feature 
around the excitation energy of 2.8 eV of the C 60 solids [12]. Coulomb interaction effects reduce 
the magnitude of the optical nonlinearity of Ceo compared with that determined using the free 
electron calculation [10], and we have discussed the possibility that the local field enhancement 
might be effective in solids. 

In this paper, we investigate nonlinear optical properties of higher fullerenes. We focus on 
the off-resonant third order susceptibility in order to estimate the magnitudes of the nonlinear 
optical responses of each isomer. The Coulomb interaction strengths are also changed in a 
reasonable range, because realistic strengths are not well known in higher fullerenes. The 
Coulomb interactions are taken into account by the parametrized Ohno potential, W(r) = 



strength at a distance r = 0, V is the strength of the long range part, and r is the mean bond 
length. Based on our results for the optical properties of C 60 and C 70 [10,11], we can assume 
V — U/2. The onsite Coulomb strength is varied within the range < U < At, t being the 
hopping integral between nearest neighbor carbon atoms. 

Figure 1 shows the absolute value of the off-resonant susceptibility, x^(0) = 0, 0, 0), 

plotted against the Coulomb interaction strength U . The different plots indicate different types 




between two electrons with distance r. Here, U is the interaction 
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of isomers. The four isomers - C 70 , D 2 -C 76 , D 3 -C 78 , and one type of C 2v -C 78 isomer [C 2v by 
Kikuchi et al's notation (Ref. 13)] - exhibit similar magnitudes of optical nonlinearities which 
are larger than those of Ceo- On the other hand, the other three isomers - another type of C 2v - 
C 78 isomer [C' 2v by Kikuchi et al's notation (Ref. 13)], D 2d -C 8 4, and D 2 -C 8 4 - show larger optical 
nonlinearities than those of the first four isomers. This is mainly due to the smaller energy 
gap of the latter isomers, even though the negative correlation between the susceptibility and 
the energy gap is not so complete through all the isomers. The decrease in the susceptibility 
between the free electron model (U = 0) and the case in which U = At is by a factor of 
approximately 0.1 for all the isomers, indicating that this is a general property of various kinds 
of higher fullerenes. The overall magnitudes of the susceptibility are around 10~ 12 esu for most 
of the Coulomb interactions considered. 

In Fig. 2, the relations between the absolute value of the off-resonant susceptibility and the 
energy gap are shown for three Coulomb interaction strengths: U = Ot, 2t, and At. Here, the 
energy gap is defined as the optical excitation energy of the lowest dipole allowed state, in other 
words, the optical gap. For each Coulomb interaction, the plots (squares, circles, or triangles) 
cluster in a bunch. When the energy gap becomes larger, the susceptibility tends to decrease. 
However, the correlation between the susceptibility and the energy gap is far from that of a 
smooth function. The correlation is merely a kind of tendency. Therefore, the decrease in the 
energy gap of higher fullerenes is one origin of the larger optical nonlinearities of the systems. 
The actual magnitudes of nonlinearities would also be influenced by the detailed electronic 
structures of isomers. 

In the calculations for C 60 reported previously, the magnitudes of the THG at the energy 
zero are approximately 1 x 10~ 12 esu in the free electron model [7], and approximately 2 x 10~ 13 
esu for U = At and V = 2t [10]. These results have been shown in Fig. 1, also. In the present 
calculations for higher fullerenes, the magnitudes are a few times larger than those of Ceo- Thus, 
the author predicts that nonlinear optical responses in higher fullerenes are generally larger than 
in Ceo- In our previous paper [10], we discussed the fact that the local field correction factor is 
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of the order of 10 for C 6 o solids. Since the distance between the surfaces of neighboring fullerene 
molecules in C 70 and C 76 solids is nearly the same as in C 60 solids, we expect that local field 
enhancement in thin films of higher fullerenes is of a magnitude similar to that in Ceo systems. 

It is of some interests to look at carbon number dependence of the magnitude of the optical 
nonlinearity of the calculated isomers in higher fullerenes. Figure 3 shows |x^(0)| as functions 
of the carbon number N for three Coulomb interaction strengths, U = Ot, 2t, and At. The solid 
lines indicate the linear fitting in the logarithmic scale: |x^(0)| ~ A ■ N a . The powers a for 
the three Coulomb interaction strengths are summarized in TABLE I. When U = Ot, the power 
a is about 5. As U increases, a decreases. It is among 4, when U ~ 2t and At. This magnitude 
of the power 4 agrees with the result of the quantum chemical calculation of higher fullerenes 
upto C 8 4 [14]. Therefore, we have shown important roles of Coulomb interactions in nonlinear 
optical response of higher fullerenes. 

Experimental measurements of optical nonlinearities in higher fullerenes whose carbon num- 
ber is larger than 70 have not been reported so much, possibly because of the difficulty in ob- 
taining samples with good quality and the difficult measurements. However, the recent report of 
the degenerate four-wave-mixing measurement of C90 in solutions [15] indicates the larger opti- 
cal nonlinearity than that in C 60 . The magnitude of is about eight times larger than in C 60 , 
and is apparently enhanced from that of the theoretical predictions: (90/60) 4 = 1.5 = 5.063. 
Therefore, further experimental as well as theoretical investigations of nonlinear optical prop- 
erties in higher fullerenes should be fascinating among scientists and technologists of the field 
of photophysics. 

In summary, we have investigated the nonlinear optical properties of higher fullerenes. The- 
oretical off-resonant third order susceptibility has been calculated using the exciton theory. 
We have found optical nonlinearities of higher fullerenes which are larger than those of C 6 o- 
The magnitude of tends to increase as the optical gap decreases in higher fullerenes. The 
nonlinearity is nearly proportional to the fourth power of the carbon number when the onsite 
Coulomb repulsion is It or At. This result, indicating important roles of Coulomb interactions, 
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agrees with quantum chemical calculations of higher fullerenes. 
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TABLE I. Coulomb interaction dependence of the power a where |x (3) (°)l ~ A-N a . 



U a 

Of 5.253 

2t 4.133 

At 3.536 
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Figure Captions 



Fig. 1. The absolute value of the off-resonant susceptibility |x^(0)| plotted against the 
Coulomb interaction strength U. The stars represent numerical results for C 60 . The closed 
and open squares represent results for C 70 and -D2-C76, respectively. The closed circles rep- 
resent results for D 3 -C 78 , and the open and crossed circles represent results for two types of 
C2D-C78. The closed and open triangles represent results for D 2 d-Cs4 and -D2-C84, respectively. 

Fig. 2. The absolute value of the off-resonant susceptibility |x®(0)| for C 60 and seven isomers 
of higher fullerenes, plotted against the energy gap (shown in units of t). The squares, circles, 
and triangles represent results for U = Ot, It, and At, respectively. The left axis is in the 
logarithmic scale. 

Fig. 3. The absolute value of the off-resonant susceptibility |x®(0)| for C 60 and seven isomers 
of higher fullerenes, plotted against the carbon number N. The squares, circles, and triangles 
represent results for U = Ot, It, and At, respectively. The left and bottom axes are in the 
logarithmic scale. The solid lines are the results of the linear fitting in the logarithmic scale: 
| x (3)| ~ A ■ N a . 
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